Fluctuations of the inverse participation ratio at the Anderson transition 
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Statistics of the inverse participation ratio (IPR) at the critical point of the localization transition 
is studied numerically for the power-law random banded matrix model. It is shown that the IPR 
distribution function is scale-invariant, with a power-law asymptotic "tail". This scale invariance 
implies that the fractal dimensions D q are non-fluctuating quantities, contrary to a recent claim 
in the literature. A recently proposed relation between D2 and the spectral compressibility \ is 
violated in the regime of strong multifractality, with \ ~~ * 1 hi the limit D2 — > 0. 

PACS numbers: 72.15.Rn, 71.30. +h, 05.45.Df, 05.40.-a 



I 

-a 
c 

o 
o 



> 

00 

o 

o 
o 
o 

-I— > 



1 

C 

O 

o 



Strong fluctuations of eigenfunctions represent one of 
the hallmarks of the Anderson metal-insulator transition. 
These fluctuations can be characterized by a set of inverse 
participation ratios (IPR) 
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In a pioneering work §, We gner found from the 
renormalization-group treatment of the cr-model in 2 + e 
dimensions that the IPR show at criticality an anomalous 
scaling with respect to the system size L, 



P q oc L- D «^-V 
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Equation (|^) should be contrasted with the behavior of 
the IPR in a good metal (where eigenfunctions are er- 
godic), P q oc L~ d ( q ~^\ and, on the other hand, in the 
insulator (localized eigenfunctions), P q oc L°. 

The scaling (||) characterized by an infinite set of crit- 
ical exponents D q implies that the critical cigenfunction 
represents a multifractal distribution Q . The notion of a 
multifractal structure was first introduced by Mandelbrot 
H and was later found relevant in a variety of physical 
contexts, such as the energy dissipating set in turbulence, 
strange attractors in chaotic dynamical systems, and the 
growth probability distribution in diffusion-limited ag- 
gregation; see for a review. During the last decade, 
multifractality of critical eigenfunctions has been a sub- 
ject of intensive numerical studies [^J. Among all the 
multifractal dimensions, D2 plays the most prominent 
role, since it determines the spatial dispersion of the dif- 
fusion coefficient at the mobility edge ||. 

In fact, to make the statement (0) precise, one should 
specify what exactly is meant by P q in its left-hand side. 
Indeed, the IPR's fluctuate from one eigenfunction (or 
one realization of disorder) to another. Should one take 
the average P q l Or, say, the most probable one? Will 
the results differ? More generally, this poses the question 
of the form of the IPR distribution function at criticality. 

In a recent Letter Parshin and Shober addressed 
this problem via numerical simulations for the 3D tight- 
binding model. Their main finding is that the fractal 



dimension D2 is not a well defined quantity, but rather 
shows universal fluctuations characterized by some distri- 
bution function ViD^) of a width of order unity. If true, 
this would force one to reconsider virtually all aspects of 
the multifractality phenomenon, such as the notion of the 
singularity spectrum f(a), the form of the eigenfunction 
correlations and of the density response at the mobility 
edge etc. In view of such a challenge to the common lore, 
the issue requires to be unambiguously clarified. 

We begin by reminding the reader of the existent an- 
alytical results concerning the IPR fluctuations. While 
the direct analytical study of the Anderson transition in 
3D is not feasible because of the lack of a small param- 
eter, statistics of energy levels and eigenfunctions in a 
metallic mesoscopic sample (dimensionless conductance 
g ^> 1) can be studied systematically in the framework of 
the supersymmetry method; see H for a review. Within 
this approach, the IPR fluctuations were studied recently 
||[To|,|g|l . In particular, the 2D geometry was considered, 
which, while not being a true Anderson transition point, 
shows many features of criticality, in view of the expo- 
nentially large value of the localization length. It was 
found that the distribution function of the IPR P q nor- 
malized to its average value (P q ) has a scale invariant 
form. In particular, the relative variance of this distribu- 
tion (characterizing its relative width) reads 



™(P q )/(P q ) 2 = Cq 2 {q-l) 2 /(3 2 g 2 



(3) 



where C ~ 1 is a numerical coefficient determined by the 
sample shape (and the boundary conditions), and /3 = 1 
or 2 for the case of unbroken (resp. broken) time re- 
versal symmetry. It is assumed here that the index q 
is not too large, q 2 <C (3ng. These findings motivated 
the conjecture || that the IPR distribution at criticality 
has in general a universal form, i.e. that the distribu- 
tion function V(P q /P q yp ) is independent of the size L in 
the limit L — > 00. Here P q yp is a typical value of the 
IPR, which can be defined e.g. as a median [|ll| of the 
distribution V(P q ). Normalization of P q by its average 
value (P q ) (rather than by the typical value P q yp ) would 
restrict generality of the statement; see the discussion 
below. Practically speaking, the conjecture of Ref. jjj 
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is that the distribution function of the IPR logarithm, 
VQxiPq) simply shifts along the x-axis with changing L. 
In contrast, the statement of Ref. Q is that the width of 
this distribution function scales proportionally to IilL. 

While the above-mentioned analytical results for the 
2D case are clearly against the statement of Q , their ap- 
plicability to a generic Anderson transition point may be 
questioned. Indeed, the 2D metal represents only an "al- 
most critical" point, and the consideration is restricted 
to the weak disorder limit g ^> 1 (weak coupling regime 
in the field-theoretical language), while all the realistic 
metal-insulator transitions (conventional Anderson tran- 
sition in 3D, quantum Hall transition etc.) take place in 
the regime of strong coupling. 

To explore the IPR fluctuations at criticality in the 
strong coupling regime, we have performed numerical 
simulations of the power-law random banded matrix 
(PRBM) ensemble. This model of the Anderson criti- 
cal point introduced in is defined as the ensemble of 
random Hermitean N x N matrices H (real for /? = 1 or 
complex for (3 = 2). The matrix elements are inde- 
pendently distributed Gaussian variables with zero mean 
(Hij) = and the variance 



(\H ij \ 2 ) = a 2 (\i-j\), 



where a(r) is given by 



a 2 (r) = 



1 sin 2 (7rr/A0 
V {ir/N) 2 



(4) 
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Here < 6 < oo is a parameter characterizing the en- 
semble, whose significance will be discussed below. The 
crucial feature of the function a(r) is its 1/r-decay for 
r 3> 6. Indeed, for r <C N Eq. (Q) reduces to 



a 2 (r) = [1 + {r/bf 



(6) 



The formula (@) is just a periodic generalization of @, 
allowing to diminish finite-size effects (an analog of peri- 
odic boundary conditions). 

In a straightforward interpretation, the model de- 
scribes a ID sample with random long-range hopping, 
the hopping amplitude decaying as 1/r with the length 
of the hop. Also, such an ensemble arises as an effective 
description in a number of physical contexts. Referring 
the reader to Refs. 12 j| for details (see also 13 - 15 ), we 
only give a brief summary of the main relevant analyti- 
cal findings. The PRBM model formulated above is crit- 
ical at arbitrary value of 6; it shows all the key features 
of the Anderson critical point, including multifractality 
of eigenfunctions and non-trivial spectral compressibility 
(to be discussed below). Perhaps, the most appealing 
property of the ensemble is the existence of the parame- 
ter b which labels the critical point: Eqs. (Q), (||) define 
a whole family of critical theories parametrized by b |L6| . 
This is in full analogy with the family of the conventional 
Anderson transition critical points parametrized by the 



spatial dimensionality 2 < d < oo. The limit b 3> 1 is 
analogous to d = 2 + e with e <C 1; it allows a system- 
atic analytical treatment (weak coupling expansion for 
the cr-model). The opposite limit 6 C 1 corresponds to 
d ^> 1, where the transition takes place in the strong 
disorder (strong coupling) regime, and is also accessi- 
ble to an analytical treatment jlTj using the method of 
ifrsf . This makes the PRBM ensemble a unique labora- 
tory for studying general features of the Anderson tran- 
sition. Criticality of the PRBM ensemble was recently 
confirmed in numerical simulations for 6=1 |l9[ . 

We have calculated the distribution function of the IPR 
in the case j3 = 1 for system sizes ranging from N = 256 
to iV = 4096 and for various values of 6 by numerically 
diagonalizing the Hamiltonian matrix defined in Eq. (^) 
using standard techniques. The statistical average is over 
a few thousand matrices in the case of large system sizes 
up to 10 5 matrices at N = 256. Specifically, we have con- 
sidered an average over wavefunctions having energies in 
a small energy interval about the band center, with a 
width of about 10% of the band width. 
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FIG. 1. Distribution of I11P2 at 6 = 1 for the system 
size N = 256, 512, 1024, 2048, and 4096. The straight 
line corresponds to the power-law asymptotics with the in- 
dex X2 ~ 4.16. 

Fig. [I] displays our result for the distribution of the IPR 
logarithm, "P^nJ^)- It is clearly seen that the distribu- 
tion function does not change its shape or width with 
increasing N. After shifting the curves along the x-axis, 
they all lie on top of each other, forming a scale-invariant 
IPR distribution. Of course, the far tail of this universal 
distribution becomes increasingly better developed with 
increasing N. From the shift of the distribution 7- > (lnP2) 
with N we find the fractal dimension D2 = 0.75 ± 0.05. 
Analogous results are obtained for other values of 6 and 
q and will be published elsewhere JTtJ . 

We conclude therefore that the distribution of IPR 
(normalized to its typical value) is indeed scale- invariant, 
in agreement with the conjecture of Ref. H and in dis- 
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agreement with Ref. H. A natural question that can be 
asked is why the authors of |?J failed to find this univer- 
sality? We speculate that, possibly, the system sizes L 
used in their numerical simulations were too small for ob- 
serving the universal form of P[ln(P 2 / P 2 VP )] m the limit 




FIG. 2. The fractal dimension D2 (squares) and the spec- 
tral compressibility \ (circles) as a function of the parameter 
b of the PRBM model. The corresponding ii > 1 and i< 1 
analytical asymptotics are shown. 

The value of the fractal dimension D 2 found from 
the scaling of the shift of the distribution with N is 
shown in Fig. || as a function of the parameter b of the 
PRBM model. The numerical results agree very well 
with the analytical asymptotics in the limits of large b, 
n = l-D 2 = 1/7T& HH and small b, D 2 = 2b 0. We 
have also calculated the spectral compressibility \ char- 
acterizing fluctuations of the number n of energy levels 
in a sufficiently large energy window 5E 1 var(n) = x{n). 
The results are also shown in Fig. |^, and are in per- 
fect agreement with the large-& asymptotics, \ — 1/27T& 
as well. A non-trivial value of the spectral com- 
pressibility < x < 1 (intermediate between \ = in a 
metal and x = 1 in an insulator) has been understood to 
be an intrinsic feature of the critical point of the Ander- 
son transition 

In a remarkable recent work [^2), Chalker, Lerner 
and Smith employed Dyson's idea of Brownian motion 
through the ensemble of Hamiltonians to link the spec- 
tral statistics with wavefunction correlations. On this 
basis, it was argued in Ref. |^3| that the following exact 
relation between \ and D 2 holds: 



X = (d - D 2 )/2d . 



(7) 



According to (]?]) , the spectral compressibility should tend 
to 1/2 in the limit D 2 — > (very sparse multifractal) , and 
not to the Poisson value \ — I. However, the numerical 
data of Fig. || show that, while being an excellent ap- 
proximation at large b (we remind that for our system 



d = 1), the relation (Q) gets increasingly stronger vio- 
lated with decreasing b. In particular, in the limit b — > 
(when D 2 — » 0) the spectral compressibility tends to the 
Poisson limit \ — > 1. The same conclusion was reached 
analytically in || for the PRBM model with broken time 
reversal invariance. Similar violation of (Q) is indicated 
by numerical data for the tight-binding model in dimen- 
sions d > 4 fl24|j . It would be interesting to see why the 
derivation ofTJ^) in [^3| fails at small b. 

Let us now comment on the necessity to distinguish be- 
tween the average value (P q ) and the typical value P q yp . 
This is related to the question of the asymptotic behav- 
ior of the distribution V(P q ) at anomalously large P q . It 
was found in the 2D case that the distribution has 



a power-law tail V{P q ) oc P q 



with x q — 2/3iTg/q 2 



(as before, g » 1 and q 2 <C f3irg assumed). We believe 
that the power-law asymptotics with some x q > is a 
generic feature of the Anderson transition point. This is 
confirmed by our numerical simulations, as illustrated in 
Fig. [I]. For not too large q the index x q is sufficiently 
large (x q > 1), so that there is no essential difference 
between (P q ) and P q yp . However, with increasing q the 
value of x q decreases. Once it drops below unity, the 
average (P q ) starts to be determined by the upper cut- 
off of the power-law "tail", determined by the system 
size. As a result, for x q < 1 the average shows a scaling 
(P q ) oc L - - ^ 1 ? -1 ) with an exponent D q different from D q 
as defined from the scaling of P q yp (see above). In this 
situation the average value (P q ) is not representative and 
is determined by rare realizations of disorder. Therefore, 
the condition x q = 1 corresponds to the point a_ of the 
singularity spectrum with /(«-) = 0. If one performs 
the ensemble averaging in the regime x q < 1, one finds 
Dq as the fractal exponent and (after the Legendre trans- 
form) the function f(a) continuing beyond the point a_ 
into the region f(a) < §. With this definition, the 
fractal exponent D q — > as q — > 00. On the other hand, 
the fractal exponent D q defined above from the scaling 
of the typical value P q yp (or, equivalently, of the whole 
distribution function) corresponds to the spectrum f(a) 
terminating at a = a_ and saturates D q — > a_ in the 
limit q — > 00. 

In the region x q > 1 (corresponding to f(a) > 0) 
the two definitions of the fractal exponents are identi- 
cal, D q = D q . This is in particular valid at q = 2 for 
the Anderson transition in 3D and for the Quantum Hall 
transition. 

As has been mentioned above, the two limits t > 1 
and 6C 1 can be studied analytically. Let us announce 
the corresponding results for the IPR statistics; details 
will be published elsewhere 0. As shown in Fig. ||, the 
"phase boundary" q c (b) separating the regimes of x q > 1 
(D q = D q ) and x q < 1 (D q > D q ) has the asymptotics 
q c = (2nb) 1 / 2 (b > 1) and q c = 2.4056 (b < 1). Notice 
that this implies D 2 = D 2 for all b. The corresponding 
power-law tail exponent x 2 is equal to irb/2 at b 3> 1 and 
to 3/2 at b <C 1. The values of x q at q < q c (for 6^1) 
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as well at q > q c are given in Fig. 
8 



2.4056 




FIG. 3. "Phase diagram" of the multifractal spectrum. 
The phase boundary q c (b) separates the regions of x q > 1 
(below) and x q < 1 (above). The dotted line separates the 
asymptotic regimes b <C 1 and b ^> 1, for which the analytical 
results have been obtained jl7| . The dashed line is a schematic 
illustration of the crossover between the two asymptotics. 

Finally, it is worth mentioning that the meaning of 
universality of the IPR distribution at the critical point 
is the same as for the conductance distribution or for the 
level statistics. Specifically, the IPR distribution does de- 
pend on the system geometry (i.e., on the shape and on 
the boundary conditions). However, for a given geometry 
it is independent of the system size and of microscopic 
details of the model, and is an attribute of the relevant 
critical theory. 

In conclusion, we have studied the IPR statistics in 
the family of the PRBM models of the Anderson tran- 
sition. Our main findings are as follows: (i) The dis- 
tribution function of the IPR (normalized to its typical 
value Pq yp ) is scale-invariant, as was conjectured in ||. 
(ii) The scaling of P* yp with the system size defines the 
fractal exponent D q , which is a non-fluctuating quan- 
tity, in contrast to |7j. (iii) The universal distribution 
V(z = P q /Pq yp ) has a power-law tail oc z^ 1 ^ Xq . At suf- 
ficiently large q one finds x q < 1, and the average value 
(P q ) becomes non- representative and scales with a dif- 
ferent exponent D q / D q . (iv) The relation (^) between 
the spectral compressibility and the fractal dimension D2 
argued to be exact in Ref. |^3| is violated in the strong- 
multifractality regime. In particular, \ — > 1 in the limit 
of a very sparse multifractal (£>2 — > 0). 

Discussions with V.E. Kravtsov, L.S. Levitov, 
D.G. Polyakov, I. Varga and I.Kh. Zharekeshev are grate- 
fully acknowledged. This work was supported by the SFB 
195 der Deutschen Forschungsgemeinschaft. 
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